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\Q Abstract 
O 

Motivated by stability questions on piecewise deterministic Markov models of bac- 
y—i terial chemotaxis, we study the long time behavior of a variant of the classic telegraph 

process having a non-constant jump rate that induces a drift towards the origin. We 
• i-H compute its invariant law and show exponential ergodicity, obtaining a quantitative 

control of the total variation distance to equilibrium at each instant of time. These re- 
sults rely on an exact description of the excursions of the process away from the origin 
and on the explicit construction of an original coalescent coupling for both velocity 
and position. Sharpness of the obtained convergence rate is discussed. 

Key words and phrases. Piecewise Deterministic Markov Process, coupling, long time 
behavior, telegraph process, chemotaxis models. 
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1 Introduction 



1.1 The model and main results 

Piecewise Deterministic Markov Process (PDMP) have been extensively studied in the last 
two decades and received renewed attention in recent years in different applied probabilistic 
models (we refer to [5] or [10] for general background). We consider the simple PDMP of 
kinetic type (Zt) t>Q = Wt)) t>0 with values in M. x { — 1, 1} and infinitesimal generator 

Lf(y,w) = wd y f{y,w) + (a + (b - a)t {yw>0} )(f(y, -w) - f(y,w)), (1) 

where b > a > are given real numbers. That is, the continuous component Yj evolves 
according to ^ = Wt and represents the position of a particle on the real line, whereas 
the component Wt represents the velocity of the particle and jumps between +1 and —1, 
with instantaneous state-dependent rate. More precisely, as long as Yj is positive the jump 
rate of the velocity is equal to b if W = +1, and it is equal to a if W = — 1; the situation is 
reversed if Yj is negative. The case a = b corresponds to the classical telegraph process in 
K x { — 1, +1} introduced by Kac [11], in which case the density of (Yj) solves the damped 
wave equation 

d 2 p d 2 p dp 
dt 2 dx 2 dt 

called the telegraph equation. The telegraph process, as well as its variants and its con- 
nections with the so-called persistent random walks have received considerable attention 
both in the physical and mathematical literature (see e.g. [9] for historical references and 
for some recent probabilistic developments). It is well known that (Yt)t>o converges when 
a = b to the standard one dimensional Brownian motion in the suitable scaling limit. 
Figure 1 shows a path of Y driven by (1) with a = 1 and 6 = 2. 

In this paper, we are interested in the long-time stability properties of the process (1) 
when b > a. One of our motivations is a better understanding of dissipation mechanisms in 
the setting of hyperbolic equations, where the telegraph process appears as the prototypical 
associated Markov process. A second motivation is to make a first step in tackling questions 
on the trend to equilibrium of velocity jump processes introduced in [6] , [7] , which model 
the interplay between intra-cellular chemoattractant response mechanisms and collective 
(macroscopic) behavior of unicellular organisms. These PDMP describe the motion of 
flagellated bacteria as a sequence of linear "runs", the directions of which randomly change 
at rates that evolve according some simple dynamics that represent internal adaptive or 
excitative responses to chemical changes in the environment (we refer the reader to [20] for 
a deeper probabilistic description) . The emergence of macroscopical drift is expected when 
the response mechanism favors longer runs in specific directions, and has been numerically 
confirmed in [6], [7]. In [20], with the aim of developing variance reduction techniques for 
the numerical simulation of these models, a so-called gradient sensing process was derived 
from them, in asymptotics where the response mechanisms of bacteria, roughly speaking, 
act infinitely fast (see Lemma 2.5 in [20] for a precise mathematical statement). The 
process (Z t ) above exactly corresponds to the gradient sensing process for the particular 
chemoattractant potential S(x) = c\x\ > in K, and constitutes a tractable toy model for 
the long-time behavior of the processes considered in [6], [7], [20] . 

When b > a a particle driven by (1) spends in principle more time moving towards the 
origin than away from it. Thus, a macroscopic attraction to the origin should take place 
in the long run, though in a consistent way with the fact that the particle has constant 
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Figure 1: Trajectory of the continuous part Y from the dynamics (1) with a = 1 and 
6 = 2. 



speed. Our main goal is to clarify this picture by determining the invariant measure \i of 
(Y, W) when 6 > a, and obtaining quantitative bounds (i.e. estimates that are explicit 
functions of the parameters a and 6) for the convergence to fi of the law of (Yt, Wt) as t 
goes to infinity. Denote by \\n — t/||tv the total variation distance between two probability 
measures rj and fj on R (recalled below at (4)). Our main result is 

Theorem 1.1 The invariant probability measure [i of (Y,W) is the product measure on 
E x { — !,+!} given by 



li{dy,dw) 



< b - a ^dy®^(8^ + S +1 )(dw). 



Moreover, denoting by fxf' w the law of Z t = (Y t ,Wt) when issued from Zq = (y,w), we 
have, for any y, y £ 1R and w, w G { — 1, +1}, 



y.w v.w 



TV 



< C(a,b)e 



r(a,b)\y\V\y\ e -\ c t 



(2) 



where 
C(a,b) 



5/2 



a + b 
/a6 + 6' 



, , . 3(6 — a) i— - 
r(a, b) = V (6 — vab) and A c 



(Vb 



We easily deduce 



Corollary 1.2 Let rj be a probability measure inMx {— 1,+1} and let ^ the law of Zt 

when the law of Zq is given by r/. Then, 

llM? - mIItv ^ °( a > h ) I e r(a ' 6)|y| (/i + rj)(dy, dw)e-^ J . 
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Figure 2: Empirical law of Y t starting at (5,-1) for t £ {2,6,10,14,18,22} with a = 1 
and 6 = 2. 

The upper bound (2) is integrable under the invariant measure fi of the full process (Y, W) 
since r(a,b) < b — a. Thus, Corollary 1.2 is significant as soon as (y,w) i-> e r{ - a ^\y\ is X]- 
integrable, ensuring in that case the convergence to equilibrium at exponential rate A c . 
Figure 2 compares the empirical law of Yt to its invariant measure for successive times 
(the shapes might be compared to those presented in [6, p. 385]). 

In spite of the simple form of the process (1), fully explicit computations on this model 
are not easy to carry out. When Yq = y > 0, as long as t < y the law of (Z s )o< s <t is equal 
to the law of the process with generator 

Hf{y, w) = wd y f(y, w) + (a + (b- a)t {w>0} ) (f(y, -w) - f(y, w), (3) 

which was computed in [9] in terms of modified Bessel functions. We have been unable 
to compute the transition laws for (1) for general time intervals. (Notice when a < b that 
long-time behavior of the process driven by (1) is completely different from that of the 
process driven by (3), which drifts to — oo.) 

The proof of the bound in Theorem 1.1 will rely on the construction of a coupling (see 
[12] and [21] for background), which classically provides a convergence rate to equilibrium 
depending on tail estimates of the coupling time. A related and popular approach to the 
long-time behavior of Markov processes is the Foster-Lyapounov-Meyn-Tweedie theory 
(see [14, 15, 19, 2]), which allows one to prove exponential ergodicity under conditions 
that are relatively easy to check. Specific applications to PDMP have been developed in 
[3, 4]. Such general results however provide convergence estimates which are not fully 
explicit, and sharpness of the bound and rates that can be deduced is hardly assessable. 

A fundamental step in proving Theorem 1.1 will be to first establish an analogous 
result for the reflected (at the origin) version of the process. A fully explicit coupling 
will be first constructed for the latter, inspired in the coalescent couplings for classic 
non- negative continuous-time Markov process, namely the M/M/l queue and reflected 
Brownian motion with negative drift. The main difficulties in our case are that we have 
to couple both position and velocity and that, contrary to those examples, we do not have 
a natural order structure. This prevents us from using the framework developed in [13] to 
construct couplings for processes that are said to be stochastically ordered. 
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We will next recall basic ideas employed to study the long time behavior of Markov 
processes via couplings, following [12]. We also introduce the "reflected version" of the 
process driven by (1) and state an analogue of Theorem 1.1 for the latter in Theorem 
1.4. The strategy of the proofs of Theorems 1.1 and 1.4 together with the structure of the 
remainder of the paper are then explained. Let us anticipate that the convergence rate A c 
in (2) will arise as the supremum of the domain of the Laplace transform of the hitting 
times of origin for the process (Yf), suggesting that this rate is sharp. In this direction, we 
also will see below that in the suitable scaling limit where (Yt) converges to the Brownian 
motion drifted to the origin, the known total variation convergence rate to equilibrium of 
the latter is recovered as the rescaled limit of the A c 's. 



1.2 Preliminaries 

c 

In all the sequel we will use the notation = meaning "equal in law to". By £(X) and V{\) 
we will respectively denote the exponential law and the Poisson law of parameter A > 0, 
whereas B(p) will stand for the Bernoulli distribution of parameter p £ (0, 1). 

Recall that the total variation distance between two probability measures rj and fj in 
a measurable space X is given by 



\V 



fj\\ Ty = inf |P(X ^ X) : X,X random variables with C(X) = r], C(X) = fj\, (4) 



where each pair of random elements (X, X) of X is simultaneously constructed in some 
probability space and is called a coupling (see [12] for alternative definitions of this distance 
and its main properties). A coupling (Ut,Ut)t>o of two stochastic processes such that 
Ut+T, = Ut+T, for any t > and an almost surely finite random time T* is called a 
coalescent coupling (T* is then called a coupling time). It follows in this case that 

c(u t )-c(u t ) < > *)■ 

A helpful notion in obtaining an effective control of the distance is stochastic domination: 

Definition 1.3 ([12]) Let S andT be two non-negative random variables with respective 
cumulative distribution functions F and G. We say that S is stochastically smaller than 
T and we write S < s to. T, if F(t) > G(t) for any t € R. 

In particular, for a couple (Ut, Ut) as above, Chernoff's inequality yields 

at\ -At 



C{Ut)-C{U t ) ^ < P(T > t) < E [e* 1 j e~ M (5) 

for any non-negative random variable T such that T* < s to. T, and any A > in the domain 
of the Laplace transform A i— > E(e AT ) of T. 

We will use these ideas to obtain the exponential convergence estimates for Z = (Y, W) 
in Theorem 1.1, and in Theorem 1.4 below for its reflected version (X, V) which we now 
introduce. The Markov process ((X t , Vt)) t>0 is defined by its infinitesimal generator: 

Af(x, v) = vd x f(x, v )+(a + (b- a)l {v>0} + (f(x, -v) - f(x, v)), (6) 

V 1 {x>0}/ 

with < a < b (the term l{x=o}(l{a;>o})~ 1 means that X is reflected at zero). The 
dynamics of the process is simple: when X is increasing (resp. decreasing), V flips to — V 
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with rate b (resp. a) and it is reflected in the origin (i.e. as soon as X = 0, V flips to 1). 
Given a path ((Y t , Wt)) t>0 driven by (1), a path of ((X t , Vt)) t>0 can be constructed taking 

X t = \Y t \, V = sgn(Y )Wo 

and defining the set of jump times of V to be 

{t > : AV t / 0} = {t > : AW t / 0} U {t > : Y t = 0}. 

Notice that since W does not jump with positive probability when Y hits the origin, 
one can also construct a path of ((Yt,Wt)) t>0 from an initial value y £ M and a path 
((X t , Vt)) t>0 driven by (6): writing gq = and (o"j) i>1 for the successive hitting times of 
the origin, we define 

{Y u W t ) = (-iy S gn(y)(X t , V t ) if t E [a u a i+1 ]. 
Let us state our results about the long time behavior of (X, V). 

Theorem 1.4 The invariant measure of (X, V) is the product measure on M+ x {—1, +1} 
given by 

v{dx,dv) = (b - a)e- ( - b - a)x dx (g) -{S-i + S+i)(dv). 

If v x,v stands for the law of (Xt, Vt) when Xq = x and Vq = v, we have, for any x, x > 
and v, v € {—1, +1}, 

(q + b)b r (a,b)(xvx) -A e i f7 x 
TV" 2a 2 ' lj 



where 



. 3(6 -a) ., AT . , A a + fe ^- (v^-v/") 2 
r ( a) 6) = ^- — ^v(6-Va6) and A c = Va6 = — 



Notice that for small times t < \x — x\/2 the total variation distance does not decrease 
exponentially fast: the distance between u t ' and v x,v is equal to 1 since the supports of 
these two probability measures are disjoint. 

Theorem 1.4 should be compared to results on two classic examples of ergodic non- 
negative continuous time Markov processes, obtained by coupling arguments that are 
briefly recalled next. Consider first Brownian motion with negative drift — c < reflected 
at the origin and which has the law £(2c) as invariant measure (see [12] for this and the 
following facts). A coupling of two of its copies (Uf)t>o and (Uf)t>o respectively starting 
from x and x consists in letting them evolve independently until they are equal for the first 
time and choosing them equal from that moment on. By non-negativity and continuity 
the coupling time T* for (Uf, Uf)t>o is stochastically smaller that the hitting time T of the 
origin for (Uf yx ) t> Q. Since for y > the hitting time T of the origin by (U^)t>o satisfies 

E y (e XT ) = exp {y{c- Vc 2 - 2A)) if A € (-oo,c 2 /2] and E y (e XT ) = +oo otherwise (see 

e.g. [17, p. 70]), taking A = c 2 /2 in (5) one gets 

\\C{Uf) - £(t/f)|| TV < e c(xV£) e- c2 ' /2 for all x,x G R+, t > 0. (8) 
This estimate can also be used to study the long time behavior of the solution of the SDE 

d£ t =dB t -c sgn(&) dt (9) 
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which has the Laplace law ce~ 2c ^dx as invariant measure. One first has to couple the 
absolute values; the first hitting time of the origin after their coupling time stochastically 
dominates the coupling time for (9). A second example is the M/M/l queue, that is the 
continuous time Markov process (Nt) t>0 taking values in N with infinitesimal generator 

Af(n) = a(f(n + 1) - f{n)) + ML {n>0} (/(n - 1) - f(n)) 

where b > a > (to ensure ergodicity). Since two independent copies of the process 
starting from n and n do not jump simultaneously and they have one unit long jumps, 
their coupling time is smaller than the hitting time of the origin T for the process starting at 
nVu. For each initial state n G N the Laplace transform of T has domain (— oo, {yb— \fa) 2 \ 
and we have E n ^e^^^ 2 ^ = (^) n 2 (see [18] for these facts) which as before yields 



C{N?)-L{N f t 



■ i \ (nVn)/2 

< ( - ) e ~(Vb-^E) 2 t for any nj ~ G Nj t > Q 



tv \a J 



We notice that in the appropriate scaling limit, the M/M/l queue is furthermore known 
to converge to the reflected Brownian motion with negative drift (see [18]). 

The construction of a coalescent coupling for the process (X, V) driven by (6) is harder 
than the previous examples since both positions and velocities must be coupled at some 
time. This will be done in two steps. In Section 3.1 we will obtain an estimate (in the 
sense of stochastic domination) for the first crossing time and position of X and X for 
a suitable coupling of the pair. At that time the velocities will be different. We will 
then construct in Section 3.2 the coalescent coupling when starting from that special 
configuration. In Section 3.3 we will obtain an explicit upper bound for the Laplace 
transform of the coalescent time, and thus the quantitative convergence bound (2). The 
required stochastic dominations will be established in terms of hitting times and lengths of 
excursions of (X, V) away from the origin. These hitting times will be previously studied in 
Section 2. We will also give therein a complete description of the excursions and compute 
thereby the invariant measure of (X, V) using a standard regeneration argument. Finally, 
Theorem 1.1 will be proved in Section 4 by transferring these results to the unreflected 
process. 

We end this section noting that A c in (2) is the right convergence rate for the process 
(1), at least in the natural diffusive asymptotics of the process. Let c > and < a N < 
b N , N G N be real numbers. We have 

Proposition 1.5 Assume that a N + b N — > oo and b N — a N — > 2c G (0, oo) as N — )■ oo. 

Let (1^ , Wf )t>o denote the process driven by (1) with coefficients a = a N ,b = b N and 
starting from a random variable Yq N ^ = ^ £ I. Then, as N — > oo, the process 

V ? * )t>0 ' V *Ow+tw)/2,/ t > 
converges in law in C([0, oo), R) to the solution of the SDE (9) with initial condition £o- 
Remark 1.6 By Theorem 1.1 and the dual representation of the total variation distance, 

\nf{t\ N),y )-nf{t {N) *))\ < c(a N , 6 JV ) e -(^-M(biviyi) exp (_ (ajv + bN ) ^-^ 2 t) 
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holds for every t > and each continuous function f : R — > [—1, +1]. Letting N — > oo and 
taking then supremum over even functions f in the previous class, we then get that 



\L{U?)-L{Uf 



I TV 



< e 



3c(xVx)/2 -c 2 t/2 



for all x, x G M+, t > 0, 



where (Uf)t>o and (Uf)t>o are Brownian motions with drift —c < reflected at the origin, 
respectively starting from x and x. Comparison with (8) suggests that the convergence rate 
A c of Theorem 1.1 cannot be substantially improved on, and that one could in principle 
improve upon the exponent r(a,b) therein (more precisely, upon the term 3(6 — a) /A). 

Proof. We will use a standard diffusion approximation argument. Omitting for a moment 
the sub and superscripts for notational simplicity, and writing j t :=W% — 2nWt(a + (b — 
a)l{Y t Wi>o})> Jt := foJsds and Yt '— Yt + K Wt for a given constant k > 0, we see by 
Dynkin's theorem that the processes Mt := Yt — Jt an d Nt '■= Y^ — t2n — JjJ 2Y s j s ds are 
local martingales with respect to the filtration generated by (Yt, Wt). Using integration by 
parts we then get that = N t — 2 J Q J s -dM s + 2nt — 2k J q W s j s ds. Thus, noting that 
j s = sgn(y s ) ((2a« — 1) + 2 x t^ Ys w s >o}(^ ~ K (b + a ))) > we see f° r K = ( a + that 



Mt = Yt 



sgn(y s — — 
o \a + b 



ds 



a + b 



Mf 



2t 
a + b 



2 / W s sgn(Y s 



(a + b) 



:ds 



are local martingales. Therefore, defining for each TV E N 



(JV) 



sgn(^))da 



W 4 



(N) 



t{a N +b N )/2 

a N + b N 



and 



(N) 



a N + 6jv 



we readily check that the processes ^ N \a)-"' and f3\"' satisfy assumptions (4.1) to (4.7) 
of Theorem 4.1 in [8, p. 354] (in the respective roles of the processes X n (t),A n (t) and 
B n (t) therein). That result ensures that C(^) converges weakly to the unique solution 
of the martingale problem with generator Gf(x) := ^f"(x) — csga(x)f'(x),f G C£°(M) 
and initial law C(£o)- D 



W 



(N) 







s(a N +b N )/2 



sgn(^)ds 



2 The invariant measure of the reflected process 

In this section we will determine the invariant measure of (X, V) . This process is clearly 
positive recurrent since, as will be shown in the sequel, the Laplace transform of the hitting 
time of (0, +1) is finite on a neighborhood of the origin, whatever the initial data are. We 
will need the following well-known results for Poisson processes. 

Proposition 2.1 ([16]) Let (Nt) t>0 be a Poisson process with intensity A > 0. Denote 
by [T n ) n>l its jump times. Then, Nt ~ "P(Xt) for any t > 0. Moreover, conditionally on 
{Nt = k}, the jump times T\, T2, . ■ . , have the same distribution than an ordered sample 
of size k from the uniform distribution on [0, t] . 
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2.1 Excursion and hitting times 

We start by computing the Laplace transforms of the length of an excursion (to be defined 
next) and of the hitting times of the origin when starting from (x, v) £ M + x { — 1, +1}. 

Definition 2.2 An excursion of (X, V) driven by (6) is a path starting at (0, +1) and 
stopped at 

S = mi{t>0 : X t = 0}. 
We denote by ip the Laplace transform of S: 

V> : AeR^(A) =E (0i+1) (e A5 ). (10) 

Notice that lim t _ 5 . >s - Vt = —1 and Vs = 1. 

Lemma 2.3 (Length of an excursion.) The domain of ifi defined in (10) is equal to 
(— oo, A c ] where 

X c = a -±±-Va-b= {V ~ b -^ )2 . (11) 
2 2 

Furthermore, if X < X c , 

= Q + 6-2A-V(a + 6-2A)^-4^ 
2o 

In particular, ip(\ c ) = \Jbja and E(o,+i)(S l ) = 2/(6 — a). 

Proof. During a time length E of law £"(6), V is equal to 1 and X grows linearly. At time 
t = E, V flips to —1 and X starts going down. Denote by Ti the second jump time of 
V. If Xt 2 = 0, then S = 2E. Otherwise, X starts a new excursion above Xt 2 which has 
the same law as S and is independent of the past. After this excursion, (X, V) is equal 
to (Xt 2 ,— !)• Once again, it reaches directly or V flips to 1 before doing so, in which 
case a new independent excursion begins. Proposition 2.1 and the strong Markov property 
ensure that, conditionally on {E = x}, the number N of embedded excursions has the law 
V(ax). We thus can decompose S as 

N 



S±2E + Y,S k , 
fe=i 

where E ~ £(&), C(N\E = x) = V(ax) and (Sk) k>1 is an i.i.d. sequence of random 
variables distributed as S and independent of the couple (E,N). As a consequence, 

V(A) =E(E^e 2XE+x ^^\E,N^ =E(e 2XE E{^(X) N \E)} 

= E ( e 2\E e aE(4>(\)-l)\ = b 

V J b + a-2X-a^(X) 

for each A in the domain of ip (which contains (— oo, 0]). This implies (since ip is a Laplace 
transform) that 

a + b - 2X - J (a + b- 2A) 2 - 4o6 

^ (A) = 2a • 

The relation is in fact valid as soon as the argument of the square root is non- negative i.e. 
as soon as A < A c with A c defined in (11). At last, E( 0)+1 \ (S) = ip'(0) = 2/(6 — a). □ 



9 



Remark 2.4 (Number of jumps in an excursion) Since each excursion is preceded 
by a jump, the number M of jumps ofV during an excursion (omitting the jump at time 
S) satisfies 

N 

M = 1 + + Mi) 
i=i 

where (Mj) i>0 is an i.i.d. sequence with the same law as M and independent of the random 
variable N such that C(N\E) = V(aE) with E ~ £{b) . By conditioning first in E,N as 
in the previous proof, one can easily derive a second degree equation and then an explicit 
expression for the Laplace transform of the number of jumps. We omit the details since 
this result will not be needed. 

Lemma 2.5 For x > 0, let S( x _i\ denote the hitting time of starting from (x, — 1). 
Then 

E(e^.- 1 ))= e -W with c{x)= b-a-V(a + b-2Xy-Aab ^ 
if A G (— oo, A c ], and +oo otherwise. 

Proof. As in the proof of Lemma 2.3, one can decompose Sy X) _i) as 

N 



k=i 



where N is a random variable with law ^(ax) independent of the i.i.d. sequence of random 
variables (<Sfc)fc>i with Laplace transform ip. Then, 



ax(ip(\)-l)+\x 



k>l 



At last, a(ip(X) — 1) + A is equal to c(A). □ 

Corollary 2.6 For any x > 0, let us denote by Sr Xi ^ the hitting time of starting from 
(x, +1). Then 

where i/j is given by (12) and c(A) by (13). 

Proof. The strong Markov property implies that S^+i) = S + Sr x _i\ where S is the 
length of an excursion independent from S( x _iy □ 

Lemma 2.7 For any x,x > 0, 

S(x+x-l) = S{x,-l) + ^(£,-1) — sto. £(a;,-l), 

where Sr x _i) and SV^-i) are independent. 

Proof. This is a straightforward consequence of the strong Markov property. □ 
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2.2 The invariant measure 

Recall that the invariant law of (X, V) is denoted by v. 

Lemma 2.8 For any bounded function / :ix { — 1,+1} — > M, we have 



where S is the first hitting time of 0. 

Proof. We will use a standard result on regenerative processes (see Asmussen [1, Chapter 
VI] for background). Let (S n ) n >i denote the lengths of the consecutive excursions away 
from 0, So := S and n := So + • • • + S n . By the strong Markov property, (0 n )neN is a 
renewal process, for each n G N the post n -process (0„+i, n +2, • • • , (Xe n +t, Ve n +t)t>o) 
is independent of (Go, • • • , @n), and is equally distributed for all n > 1. This means that 
(X t ,Vt)t>o is a regenerative process with regeneration points (0 n )neN and cycle length 
corresponding to the length of an excursion. The result is immediate from [1, Theorem 
1.2, Chapter VI] and Lemma 2.3. □ 

Lemma 2.9 Define, for a non negative function g : { — 1, +1} —> K and A G K, 

rS 



[ e xx °g(V s )ds. 
Jo 



Then, conditionally on (Xq, Vq) = (0, +1), we have 

r E N f s(i) n 

F= (g(l) +g(-l)) / e Xy dy + y2 e KEU^ N) +x s ) g( yP)ds, where 
Jo i=1 Jo 

• (^i)i>o ^ S an se Q uence °f independent uniformly distributed random variables on 
[0,1], and for each n > 1 (Un n \, Ur 2>n y ■ ■ ■ , Ur n>n \) is the re-ordered sampling of 
{U u U 2 ,...,U n ); 

• (xjp ,Vt) 0<t<S (i) is a sequence of independent excursions; 

• E ~ £ (b), C(N\E = x) = V{ax), and the pair (E,N) is independent of the all the 
previous random variables. 

Proof. The argument has been already been given in the first part of the proof of 
Lemma 2.3. We just notice that the N independent embedded excursions therein occur 
at the heights (EUr N m, EUr N _i jN \, . . . ,EUnj,n) (see Proposition 2.1). □ 

We now are ready to compute v which is the first point in Theorem 1.4. Since 

E (i) 5 eA(B ^ )+X " l))5 (^ (i) ) ds | jE ' Ar ' (C/4) ^ 1 ) = eW(l,iV)IE (o,+i)( F )' 
we get that 

E (0 ,^(F) = e( feffl + _ i) + Em ^ (F) f^ M V 



V i=i 



n 



We have, for any x > and k E N, 



k / k \ 

^ E (e A ^>) = E [ ^ e XxU > J 

i=l ' \i=l / 



iV 



and then 

E ( ^2 e XEU V> N) 
\i=i 

As a conclusion, we have 



Ax 



1). 



= X V^A " J 

g(l) + g(-l) a 

" 6 - A +6TA E <°.+i)( n 
which provides the expression of Em +i) (^) for any A < 6 — a: 

<7(1) + <?(-!) 



E ( o,+i)(^) 



6 — a — A 



On the other hand, Lemma 2.3 ensures that E( 0) +i)(<S l ) = 2/(6 — a) from where we get, for 
any A < b — a, that 



/ e Xx g(v)v(dx, dv) 

iRx{-l,+l} 



b-o g(l) + g(-l) 
b—a—X 2 



In other words, the invariant measure of (X, V) is v = £(b — a) <g> (l/2)(<5_i + c) + i} 



3 The coalescent time for the reflected process 
3.1 The crossing time 

We will first construct a coupling (X,V,X,V) starting at (x,v,x,v) until a time T c = 
T c (x,v,x,v) called crossing time, at which Xt c = Xt c - In doing so, we will also stochas- 
tically control T c and Xt c - The coupling will consist in making the two velocities equal 
as longer as possible. Assume without loss of generality that x < x. Plainly, if V and V 
are different, we let the two processes evolve independently until one of them performs a 
jump or until X — X hits 0. At that time, if X ^ X, the two velocities are equal and 
we set them equal until X hits the origin. During this period the paths of X and X are 
parallel and, at the hitting time of the origin, V and V are once again different. We then 
iterate this procedure until T c . Notice that X is smaller than X on [0, T c ). We now make 
the construction with full details. 



3.1.1 The main initial configuration 

Assume first that (Xq, Vq, Xq, Vq) = (x, — 1, 0, +1). The coupling works as follows: with 
rate a (resp. 6), V (resp. V) flips to +1 (resp. —1) and if none of these two events occurs 
before time x/2, then 

X x /2 = X x j 2 and V x / 2 = -1 = —%/%■ 
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If a jump occurs at time t\ < x/2, then (X Tl , V Tl , X T1 , V n ) = (x — t\,U,t\,U), where 
U = —1 with probability b/(a + b) (V jumps before V) and U = 1 with probability 
a/(a + b). Then, V and V are chosen equal until X hits z.e. during a time Sr T1 m and 

(^ri+5 (Tlia) ,Ki+S (n ^ )J ^ri+% 1]C/) ,Ki+% l!(7) ) = (x- 2ti, -1,0, +1). 

Notice that S7 T1)+1 ) = 5 + S'fo i) where 5 is the length of an excursion independent 
of Sfa \y As a conclusion, if a jump occurs at time T\ < x/2, then the full process 
(X, V, X, V) is equal to (x-2ti, -1, 0, +1) at time n+S^^+BS where B ~ B(a/(o+6)). 
One has to iterate this procedure until X — X hits 0. 

Consider now a Poisson process (N(t)) t>0 with intensity a + b. We denote by (T n ) n>0 
its jump times (with To = 0) and define (Ti) i:>1 by T{ = Ti — Tj_i for i > 1. The number of 
return times at for X before T c is distributed as N(x/2) and the length of the periods 
when (X, V) and (X, V) are independent are given by t\, T2,. . • ,Tjv(a:/2) an d x/2 — T N r x /2\- 
Then, 

iV(a/2) 

T c (x,-1,0,+1) £ J] + + S^W) + x/2-T N{x/2) , 

i=l 

where the law of (5^)^ is the one of the length of an excursion, (S , ( Tij _ 1 )) i>1 are the 
hitting times of starting from (t*)"^, (Bi) {>1 have the law B(a/(a + 6)) and all these 
random variables are independent. Since T N ( x / 2 ) = t\ + t 2 + ■ ■ ■ + tj^i x i 2 ), Lemma 2.7 
ensures that 

JV(a;/2) 

r c (x,-i > o I +i)4 aj /2 + 5 PW( . /ajf _ 1) + X; (^w) 

i=l 

< sto . x/2 + S^/a,-!) + E(x) (14) 

where 

JV(z/2) 

S(x):4 ^ (^«). 
i=l 

Notice that E(ii + v) = + S(u), where X(u) and X(t>) are independent, and that 
distributes as the sum of the lengths of N independent excursions, with N ~ V(au/2). 

3.1.2 Other configurations 

We next construct the paths until T c {x,v,x,v) and control this time irrespective of the 
initial velocities. Without loss of generality, we can assume that x > x. We just have to 
construct the paths until (X, V, X, V) reaches a state (it, —1, 0, +1), and then make use of 
the previous section. 

Assume firstly that v = v = U G {—1, +1}- We have to construct a trajectory of (X, V) 
until S(x,U) the hitting time of 0. Define for any t £ [0, Sr x m), Vt = Vt, X t = X t — x + x, 
Vs { . v) = — 1 and Xs^ V) = x — x. Using Lemma 2.7 and (14), one has 

T c (x, U, x, U) = %,c/) + T c (x - x, -1, 0, +1) 

= St{ U=+ i} + % -i) + T c (x - x, -1, 0, +1) 

<sto. St{ U=+ iy H h ^((x+ij/a.-i) + ^{x - x). 
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Assume now that v = 1 = —v. The processes (X, V) and (X, V) are chosen indepen- 
dent until the first jump time. This is equal to E = (E\ A x) A E 2 , where E\ ~ £(a), 
E 2 ~ f (6) and Vg, X B , V") = (x + E, U, x - E, U) with U G {-1, +1}. In particular, 
one has 

x — x 



X E + X E x + x , - X E 

and 



2 2 2 2 

Since for any y,y>0, T c (y, -1, y, -1) < sto . T c (y, 1, y, 1) this ensures that 

T c (x,+l,x,-l) < sto .E + T c (x + E,+l,x-E,+l) 

x — X 



< sto . E + S + 



+ E + S^ x+x y 2 -i) + S(x - x + 2E) 



x — X 



If t> = -1 



< sto . IE + E(2E) + S + + ^(^/a _ 1} + £(x - x). 

-v, we proceed as in the previous case. With the same notations, 

X E + X E x + x \X E - X E \ x - x 
~ = — ^ — and J < — h E. 



We then get the same upper bound as before. As a conclusion we have established the 
following upper bound for T c : 

Lemma 3.1 For any x > x and v, v G { — 1, +1}, 

T c (x, v, x, v) < sto . 2E + V(2E) + S + S {{x+x)/2 _ 1} + + S(x - x), 

where E = F A x mi/t F is an exponential variable with parameter a + b and is t/ie 
sum 0/ i/ie lengths of N independent excursions where N ~ V(au/2). Moreover, 

X Tc = X Tc < + £ and V Tc = -V Tc . 



3.2 A simple way to stick the paths 

We now assume that (Xo,Vo) = (x, 1) and (Xq,Vo) = (x, — 1) and construct two paths 
which are equal after a coalescent time T cc {x). The idea is to use the same exponential 
clocks for both paths but in a different order. We explain the generic step of this construc- 
tion considering R and Q two given independent random variables with respective laws 
£ (a) and £ (b). There are two possible situations: 



Case 1: R < x. In this case, defining T = R + Q 

'+1 ift€[0,Q), 
V, =1-1 if t € [Q,T), and Vt = < 
+1 if t = T 



'-1 iftG[0,i?), 
+1 j£t€[R,T), 
-1 if t = T, 



one has Xt = x + Q — R = Xt and Vr = 1 = — Vr- 
• Case 2: R > x. In this case, defining T = x + Q, 



V t 



+1 ift€[0,Q), 
-1 if te[Q,T), 



and Vt = < 



'-1 ift€[0,x), 
+1 i£t€[x,T), 
-1 if t = T, 
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Figure 3: Two paths starting at x = 3 with different velocities until they stick together. 



one has Xt = Q = Xt and Vr = Vr = — 1. In this case (X, V) and (X,V) are 
coupled at time T. 

We now construct the paths. We take an i.i.d. sequence of independent pairs of exponential 
variables (R n ,Q n ) with R n ~ £{a) and Q n ~ £(b), and inductively define tq = and 
r n+ i = r n + T n , with T n defined from (R n ,Q n ) as above until Case 2 occurs. At each 
iteration, 

• if X does not hit the origin in the interval [r n ,r n +i] (Case 1) we set 

X Tn+1 = X Tn+1 and V Tn+1 = 1 = -V Tn+1 ; 

• if X hits the origin in the interval [T„,r n+ i] (Case 2) we set 

x r n +i = X Tn+1 , V Tn+1 = V Tn+1 = -1 and T cc (x) := r n+1 . 

By construction, Xt > Xt for any t > and the coupling time T cc (a;) is smaller than 
the hitting time of the origin time of X (see Figure 3). As a conclusion we have shown 
the following result. 

Lemma 3.2 There exists a coupling of (X, V) and (X, V) starting respectively from (x, 1) 
and (x, —1) such that the coalescent time T cc {x) is (stochastically) smaller than Si x ^\\ and 

X S (X , +1) = X S (X<+1) = and V S(x+1) = V S(X!+1) = 1. 
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3.3 The Laplace transform of the coupling time 

We now gather the previous estimates to control the Laplace transform of the coupling 
time of the two paths starting respectively from (x,v) and (x,v): 

Proposition 3.3 For any x > x > 0, any v,v £ {— 1,+1}, there exists a coalescent 
coupling such that the coupling time T(x, v, x, v) is stochastically smaller than a random 
variable 

T(x, x) = F + £(2F) + + S + S + S (X) _i) + + £(x - x), (15) 

where F ~ £ (a + b), S and S are excursion lengths and all the random variables are 
independent. Furthermore, for any A £ [0,A C ] ; 

4^' 5) ) < 2a+b ^l <X) - («w + ^+ x ^m-r 

At last, a realization of T(x, x) is the first hitting time at of X after T c (x, v, x, v) + 
T C c(X Tc(x ^), and then X T{x ^ = X T{x ^ = and V T{x ^ x) = V T(x ^ x) = 1 hold. 

Proof. From the previous sections, one can construct a coalescent coupling with a coa- 
lescent time T such that 

T(x, v, x, v) = T c (x, v, x, v) + T cc (X Tc<yX ^^^))- 



Thanks to Lemmas 2.7, 3.1 and 3.2, we get that 

x 

i\ -l — 

2 



x — X 

T(x, v, x, v) < sto . 2E + T,(2E) + S + S {(x+s .) /2 -i) H 7. h S^-a-^+E.+l) + - 



< sto . 2£ + £(2£) + +S + S + S (SB _!) + + £(x - x), 

where S is an independent copy of S and all the random variables of the right hand side 
are independent. Recall that E is equal to F A x where F is a random variable variable of 
law £{a + b). In particular, 2E < s to. F -\- x and then 

T(x,v,x,v) < st0 . T(x,x), 

where T(x, x) is given by (15). At last, for any A < A c , one has 

A + a(ip(X) - 1) + c(A) < A c + a(ip(X c ) - 1) + c(A c ) = 6- a<a + 6. 

This ensures that for any A < A c , 

a + b 

2a + b- A- a^(A) - c(A)' 
Using the independence of the random variables provides the desired upper bound. □ 

Corollary 3.4 /n particular, if x < x, 
where X c is given in (11) and 

r(a,&)= 3(6 7 Q) V(b- \/ab). 



E ^ e AF+AS(2F)+AS (F ,_ 1) ^ _ E ^(A+a^(A)-l)+c(A))F^ 
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Proof. Let us choose A = A r . Since 



X c = iVb ^, <X c ) = b -^ and ^{K) = A 



we get that 



xc{Xc) + + ^a(V(A c ) - 1) = ^ - ^ ( 3x(Vb + y/B) + x(Vb - 3^) 



If x is fixed, then the right hand side is a linear function of x G [0, x] and it is bounded 
above by the maximum of its values at x £ {0, x}. In other words, 

xc(A) + A + ^^a(^(A) - 1) < ^(6 - a) V (b - Vab). 

which concludes the proof. □ 
Using inequality (5), the end of the proof of Theorem 1.4 is now obvious. 



4 The unreflected process 

We finally sketch the proof of Theorem 1.1. The invariant measure is obtained by a 
similar regeneration argument as the one in Lemma 2.8, using the obvious relation between 
excursions away from (0, +1) of the reflected and unreflected processes, and Lemma 2.9. 
The sketch of the proof of the bound (2) is the following: 

• construct a coupling (X, V, X, V) starting from (x, v, x, v) until time T(x, x), where 

x = \y\, v = sgn(y)w, x = \y\, and v = sgn(y)w, 
and notice that Xj,^ x = X^, ~\ = (see Proposition 3.3); 

• construct (Y,W,Y ,W) on [0, T(x, x)\ from (X, V, X,W) and (y,w,y,w) (see Sec- 
tion 1.2). Notice that Y T(x ^ = Y T{x ^ = 0, but in general W T{x£) = -W T(x ^; 

• wait for the first jump time E ~ £(2b) of (Y, W, Y,W) (as the minimum of two 
independent random variables of law £ (6)); 

• construct a coalescent coupling (1", W, Y, W) starting from (E, w, —E, w) with a cou- 
pling time smaller than the hitting time of the origin when starting at (E, +1). 

We just give the details of the last point, the other ones being clear. The construction 
is similar to the one of T cc (x) for the reflected process. Assume that y = —y > and 
w = w = +1 and consider two independent random variables (R, Q) with respective laws 
£(a) and £(b). Then we may have: 



Case 1: R < y. In this case, defining T = R + Q, 

Wt -- 



+1 i£t€[0,R), 
-1 if t€[R,T), and W t 
+1 if t = T, 



'+1 ift€[0,Q), 
-1 ]£t€[Q,T), 
+1 if t = T, 



one has Yp = y + Q — R = —Yp and Wp = Wp = 1. 
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Case 2: R> y. In this case, denning T = y + Q, 



W t 



f+i if t e [0,i/), 

+1 ift€[y,r), and W t 
-1 ift = T, 



+1 ifi€[0,Q), 
-1 i£t€[Q,T), 



on has It = Y T and Wt = W T = -1. At that time, (Y, W) and (Y, W) are coupled. 

The algorithm to construct the paths (Y, W, Y, W) consists in repeating the above con- 
struction until Case 2 occurs for the first time. This will happen before Y reaches the 
origin. From this scheme and previous work on the process (X,V), the coupling time 
S(y, w, y, w) satisfies 

S(y,w,y,w) <sto. T(\y\, \y\) + E + 5 (E)+1) = T(\y\, \y\) + E + S + S (E _iy 

As a conclusion, E(e A5 ^' w '^) < E (e xT ^'^Ai/j(X) 2b _x- c (\) • In particular, 

E ( e X c S(y,w,y,w)\ < f b\ 5/2 a + b &T ^ X 

V / - \a) ^b + b 

3(b-a) 



where r(o, b) = 1 4 a ^ V (6 — yab). Using (5) ends the proof. 
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